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A Model of Polarized X-ray Emission from Twinkling 
Synchrotron Supernova Shells 
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ABSTRACT 

Synchrotron X-ray emission components were recently detected in many young super- 
nova remnants (SNRs). There is even an emerging class - SN1006, RXJ1713. 72-3946, 
Vela Jr, and others - that is dominated by non-thermal emission in X-rays, also prob- 
ably of synchrotron origin. Such emission results from electrons/positrons accelerated 
well above TeV energies in the spectral cut-off regime. In the case of diffusive shock ac- 
celeration, which is the most promising acceleration mechanism in SNRs, very strong 
magnetic fluctuations with amplitudes well above the mean magnetic field must be 
present. Starting from such a fluctuating field, we have simulated images of polarized 
X-ray emission of SNR shells and show that these are highly clumpy with high po- 
larizations up to 50%. Another distinct characteristic of this emission is the strong 
intermittency, resulting from the fluctuating field amplifications. The details of this 
"twinkling" polarized X-ray emission of SNRs depend strongly on the magnetic-field 
fluctuation spectra, providing a potentially sensitive diagnostic tool. We demonstrate 
that the predicted characteristics can be studied with instruments that are currently 
being considered. These can give unique information on magnetic-field characteristics 
and high-energy particle acceleration in SNRs. 

Key words: radiation mechanisms: non-thermal — polarization — X-rays: ISM — 
(ISM:) supernova remnants — shock waves. 



1 INTRODUCTION 

Electrons and positrons accelerated to TeV energies by dif- 
fusive shock acceleration (DSA) in SNR shells will effi- 
cient ly radiate in X-rays in the associated magnetic fields 
(e.g., [Reynolds fc Chevalieilll98l|V A recent review on 'SNRs 
at high energy' is given by Reynolds! (|2008D, In some 
sources (e.g. SN1006, RXJ1713. 72-3946 and Vela Jr), the 
synchrotron component is dominating the X-ray emission, 
whereas in others such as Cas A it is not easy to distinguish 
the synchrotron component from the bremsstrahlung emis- 
sion. Mapping of the polarized X-ray emission from SNRs 
would allow to separate out and study the synchrotron com- 
ponents. 

With the high-resolution imaging capability of Chan- 
dra, likely synchrotron structures a re already seen in the 
X-ray images of various SNRs (e.g., IVink fc Lamind 120031 ; 
iBamba et al.ll2006l ; iPatnaude fc Fesenl 120091 )." The observed 
non-thermal emission is concentrated in very thin (arc- 
seconds width) filaments and clumps and has typically a 
rather steep spectrum with an exponential roll-off. In ad- 
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dition, lUchivama et alj (|2007l ) reported variability of such 
X-ray hot spots in the shell of SNR RXJ1713. 72-3946 on 
about a one-year timescale. The thin filamentary structures 
can be naturally explained in the DSA scenario: options 
are 1) a narrow spatial extend of the TeV-regime electron 
population caused by efficient electron cooling due to syn- 
chrotron energy losses in the vicinity of the SNR shock with 
strong magnetic-field amplification (e.g., IVink fc Lamind 
l2003t IBamba et alj|200l iVinkl l200ct) and 2) the observed 



narrow filaments are limited by magnetic field damping and 
not by the ener gy losses of the radiating electrons (e.g., 
iPohl et al.ll2005l h 

Polarized X-ray emission - from any source - was 
observed so far only in very few cases. Observations of 
the Crab Nebula with X-ray polarimeters aboard OSO-8 
(|Novick et all Il972l ; IWeisskopf et all Il976t ) revealed a po- 
larized flux of about 15% at few keV energi es (also de- 
tected with the I BIS detector on IN TEGRAL bv lForot et~ai1 
^OOct l). Recentlv. lGotz et all (|2009h reported variable polar- 
ized emission at 200-800 keV from GRB 041 21 9A with IBIS. 
Very little else can be reported thus far. 

Efficient DSA of protons and electrons in supernova 
shells requires turbulent magnetic fields, with energy den- 
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sities that are a substantial fraction of the shock ram 



press ur e (e.g. , iBlandford fc Eichlerlll987l; iMalkov fc Druryl 
200 ll: Iffillai 120051: iBelll |2004 : lAmato fc Blasil 12009 : 



Vladimirov et all 20061 ) . Both regular and stochastic mag- 



netic fields determine the spectra and maps of synchrotron 
radiation of high-energy electrons and positrons from SNRs. 

A model of non-thermal radio emission from SNRs, ac- 
counting for the orientation of the r egular ambient mag - 
netic field, was presented recently by IPetruk et al.l (|2009h . 
These authors synthesized radio maps of SNRs, making var- 
ious assumptions on the dependence of the electron injec- 
tion efficiency on the shock obliquity. Their method uses 
the azimuthal profile of the radio surface brightness as a 
probe of the orientation of the ambient magnetic field. The 
effect of random magnetic fields in super nova shells on ra- 
dio sy nchrotron emission was addressed bv lStroman fc Pohll 
(|2009l ). They discussed the emission and transport of po- 
larized radio-band synchrotron radiation near the forward 
shocks of young shell-type supernova remnants with a strong 
amplification of the turbulent magnetic field. Modeling the 
magnetic turbulence was done as a superposition of waves at 
a particular moment in time; no time evolution was consid- 
ered. They found that isotropic strong turbulence produces 
only weakly polarized radio emission even in the absence of 
internal Faraday rotation. If anisotropy is imposed on the 
magnetic-field structure, then the degree of polarization can 
be significantly increased, if the internal Faraday rotation is 
inefficient. 

It has long been known that random directions of mag- 
netic fields in addition to Faraday rotation may strongly 
reduce the a verage polariza t ion of synchrotron emission 
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IStroman fc Pohlll2009l l. This explains the relatively low po 
larization frequently observed for radio synchrotron sources. 
However, as we will show below, the turbulent magnetic 
fields that reduce the average polarization can result in 
highly polarized patchy structures potentially observable in 
hig h resolution imag es at X-rays. 

lRevnold"sl (|l998l ) simulated X-ray synchrotron images 
assuming a regular magnetic field and distributions of ultra- 
relativistic electrons accelerated by a forward shock using 
age-limited and loss-limited parameterizations. 

The effect of turbulent magnetic fluctuations (includ- 
ing field magnitude fluctuations) on syn chrotron X-ray emis - 
sion images was recently addressed bv lBvkov et al. I (|2008h . 
A system of finite size filled with a random magnetic field 
was modeled and used to construct synchrotron emission 
maps of a source with kinetically simulated distributions of 
ultra-relativistic electrons. The random field was composed 
of a superposition of magnetic fluctuations (transverse plane 
waves propagating with some phase velocity) with random 
phases and a given spectrum of amplitudes. Accounting for 
the field magnitude fluctuations was especially important in 
view of the dependence of the emissivity on the local mag- 
netic field (further addressed below). A particulary strong 
dependence occurs in the cut-off regime of the synchrotron 
spectrum (also further addressed below) jBykov et alj (|2008l ) 
found that non-steady structures (dots, clumps, and fila- 
ments) typically arise, in which the magnetic field reaches 
exceptionally high values. These magnetic-field concentra- 
tions dominate the synchrotron maps, with an evolving, in- 
termittent, and clumpy appearance. The modeling showed 



that the overall efficiency of synchrotron radiation from the 
cut-off regime of the electron spectrum can be strongly en- 
hanced in a turbulent field with some \J (B 2 ), compared to 
emission from a uniform field of the same magnitude sj (B 2 ), 
but of just a random direction. Strong temporal variations 
of the brightness of small structures were found, with time 
scales much shorter than variations in the underlying par- 
ticle distribution. The variability time scale depends on the 
phase velocity and the spectrum of magnetic fluctuations. 
The simulated structures indeed resemble the 'twinkling' 
structures that are observed in X-ray images of some su- 
pernova remnants. 

The same electrons that are producing X-ray syn- 
chrotron emission will emit TeV photons by inverse- 
Compton scattering. Both processes are of fundamental im- 
portance for our understanding of high-energy particle ac- 
celeration and the distinction between leptonic and nucle- 
onic contributions to the obs erved gamma-ray emission (e.g., 
lAharonian et al.ll2007l , 120091 ). Gamma-ray images of a SNR 
with efficient DSA in differen t circu mstellar environments 
were constructed bv lLee et al.l (|2008l ). 

In this paper we expand upon the work of lBvkov et al.l 
(|2008l ). modeling the specific features of the polarized syn- 
chrotron emission arising from the stochastic nature of mag- 
netic fields of young SNR shells. In [J2]we describe the simu- 
lation setup that includes the kinetic model of a TeV regime 
electron distribution and a simulation of a random magnetic 
field with different fluctuation spectra. In [J3]we present the 
resulting polarized synchrotron emission maps for different 
X-ray energies and different magnetic fluctuation spectra. 
In [J4] we discuss the observational perspective. 



2 THE MODEL 

In order to construct maps of polarized synchrotron emis- 
sion from SNR shells, it is convenient to use the local densi- 
ties of the Stokes parameters. Because of the additive prop- 
erty of the Stokes parameters I, Q, U, V for incoherent pho- 
tons, we can integrate these over the line of sight weighted 
with the distribution function of radiating particles. The 
degree of polarization is determined in a standard way as 

n = VQ 2 + u 2 + v 2 /i. 

The synchrotron emission is characterized by a coher- 
ence length I / that is of the order of a M eV electron gyro- 
radius (see e.g. lRvbicki fc Lightmanlll979l ). In the simulation 
we only consider the effects of magnetic fluctuations having 
scales that are much larger than // ~ m e c 2 je\J (B 2 ). This 
is because in the nonlinear DSA modeling of non-relativistic 
SNR shocks the magnetic fluctuation spectra are expected 
to fall down steeply at spatial scales belo w the gyro-radius 
of a G eV proton (see for instance Fig. 3 in I Vladimirov et al.l 
I2006T ). That means that the fluctuation wavenumbers k sat- 
isfy fc • If <C 1. Therefore, neglecting the magnetic fluctu- 
ations of the scale less or c omparable to If, we apply the 
standard formulae (see e.g. iGinzburg fc Svrovatskiil Il965l ) 
for the synchrotron power of a single particle of Lorentz fac- 
tor 7 ^> 1 in the simulated random magnetic field composed 
of the long-wavelength MHD fluctuations. Then we integrate 
this power over the line of sight through the system filled 
with random field fluctuations. 

The modeling of the polarized synchrotron emission 
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Figure 1. Geometry of the simulated supernova shell. The left 
panel shows half of the shell quarter with the boxes in which 
the random magnetic field was simulated. The local densities of 
the Stokes parameters were then integrated over the line of sight 
along the axis — oo < z < oo (i.e. the two half-quarters of the 
sphere). The right panel is the resulting projection that is shown 
in the simulated maps. 



from relativistic shocks of GRBs, pulsar wind nebulae and 
AGNs objects would likely require strong small scale mag- 
netic fluctuations of wavenumbers k ■ If ~ 1 and will be 
discussed elsewhere. The first particle-in-cell simulations of 
relativistic shock s in unmagnetize d electron-positron pair 
plasmas (see e.g. ISpitkov skv 2008]) have demonstrated the 
feasibility of self-consistent modeling of pair acceleration to 
energies above 100 times that of the thermal energy. The 
simulated nonthermal particles were carrying about 10% of 
the downstream thermal energy, promising potential appli- 
cations to the modeling of polarized synchrotron emission 
from GRBs, blazars and pulsar wind nebulae. 

We start with the spectral flux densities p^\0, r y) and 
pffl(Q, 7) with two principal directions of polarization ra- 
diated by a particle w i th Lo rentz factor 7, as given by 
iGinzburg fc Svrovatskiil (|l965l ) [their Eqs.(2.20)]. Here 6 is 
the angle between the local magnetic field B(r, t) and the 
direction to the observer. In the case of a random magnetic 
field it is convenient to use the local spectral densities of the 
Stokes parameters expressed through piP and p^: 
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where the angle \ is between the major axis of the po- 
larization ellipse and a coordinate in the plane perpen- 
dicular to the observer direction, and tan f3 is determined 
by the ratio of the minor an d major axes of the ellipse 
jGinzburg fc Syrovatskiill 19651 ). 



2.1 Random magnetic field and particle 
distribution 

Synchrotron X-ray emission is radiated by 10 TeV 
regime electrons since the magnetic field amplitude in 
SNR shells is typically below a mG. Efficient DSA 



of high-energy particles requires a substantial ampli- 
fication of magnetic field fl uctua ti ons in the vicinity 
of the shock; see e.g. iBelll Jl978h : iBlandford fc Eichlerl 
l|l987l ); iMaikov fc Drurvl (|200lJ ). Magnetic-field amplifica- 



tion mechanisms due to cosmic- ray instabi li ties in nonlinear 
DSA w e re proposed recent ly by [ Belli d20o3) ; I Amato fc Blasil 



IVladimirov et al.l (|2008f k IZirakashvili fc Ptuski 



J2006I); IVladimirov etaD (|2006|); iPelletier et al 
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The models predict amplified magnetic-field amplitudes well 
above the interstellar field in the far upstream of the shock. 
These current models are suited to estimate the amplitudes 
and spectra of amplified magnetic fluctuations averaged over 
some macroscopic spatial and temporal scales. The averaged 
magnetic-field spectra available from the models are appro- 
priate to model energetic particle spectra, but do not allow 
to simulate the synchrotron images of SNR shells and to 
judge about their temporal evolution. 

In reality the distribution of the emitting electrons is 
a random function of position and particle energy because 
of the stochastic nature of both the electromagnetic fields 
and the particle dynamics. However, no self-consistent treat- 
ment of such a particle distribution in strong magnetic tur- 
bulence is available. Rigorous modeling of the magnetic-field 
structure and evolution should invoke simultaneously fully 
nonlinear PIC-type simulations of the collisionless shock, 
supersonic flow, and the effect of the high-energy parti- 
cles. A microscopic selfconsistent description of magnetic- 
field fluctuations that are strongly coupled with electric cur- 
rents of accelerated particles is not feasible yet for non- 
relativistic shock simu lations in SNRs (see appendix in 
IVladimirov et al]l2008l . for a discussion). 

Therefore, to model the synchrotron SNR images we 
simulated local statistically stationary random magnetic 
fields of giv e n spe ctra using the technique described in 
iBvkov et al.l (|2008f ). The statistically isotropic and homo- 
geneous random field was constructed as a sum over a large 
number of plane waves with wave vector, polarization, and 
phase chosen randomly. In the simulation presented below 
we assume a plane wave frequency w„(k n ) = u p h • k n pa- 
rameterized with a phase velocity u p h. The spectral energy 
density of the magnetic-field fluctuations of wavenumber k 
is described as W(k) oc k~ s , where 5 is the spectral index. 



(B 2 ) = 



dkW(k), 



The average square magnetic field {B 2 }, the spectral index 
8, and the wavenumber range (fc m i n , fc max ) are the input pa- 
rameters of the model. The spectral index 5 in the standard 
DSA scenario is expected to be in the range 1 ^ <5 ^ 2. 

Then the local spectral emissivity of polarized syn- 
chrotron emission was determined at various times using 
a calculated model distribution of electrons (or positrons). 
The kinetic model used to si mulate the electron distribution 
was described in detail by IBvkov et al.l (|2000l ). The spa- 
tially inhomogeneous electron distribution function is cal- 
culated from the kinetic equation for electrons at a SNR 
shock that uses piece-wise parametrization of the particle 
diffusion coefficient to account for both Fermi I and II types 
accelerations and that is consistent with the magnetic fluc- 
tuation spectrum. The model assumes a diffusion coefficient 
k(p) oc p a , where a — 1 for TeV-regime electrons (the Bohm 
type diffusion regime) and it has a flatter energy dependence 
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at MeV regime energies. The synchrotron losses of 10 TeV 
regime electrons in magnetic fields of \/ < B 2 > > 10~ J G 
are faster than the inverse Compton losses that are domi- 
nated by CMB photon scattering. 

The scale sizes of the particle distributions of the shock 
upstream (both electrons and protons) for DSA is A" w 
K(p)/u s h ~3x 10 17 • u s h8 • B~q ■ Er e v cm, where the r.m.s. 
magnetic field B^o is in /iG units, u s hs is the shock ve- 
locity in units of 1,000 km s" 1 , and Erev is the electron 
energy in TeV units. In the shock upstream the width of the 
layer where the highest energy electrons are stopped due to 
synchrotron losses is about A™. In the shock downstream 
the width of the ultra-relativistic electron/positron cooling 
layer A d is about A d ~ 6 x 10 21 • u shg • B~<i ■ E^ v cm. Both 
widths A u ' d of the electron regions emitting X-rays are rel- 
atively narrow, typically below 0.3 pc for B^g > 30 and 
EtcV 3> 1. Therefore, for large enough SNR shells of radii 
iW » A"' d the one dimensional approximation for the 
determination of the distribution function is well justified. 
It is also important that the wavelengths of magnetic fluc- 
tuations in the SNR shell are of the order of the gyroradii of 
the relativistic protons in DSA models, and therefore that 
these are below A"' d justifying the use of a homogeneous 
r.m.s. field in the losses term of the kinetic equation. We 
numerically calculated the electron distribution function in 
the vicinity of the SNR forward shock. The results were then 
used in simulations of the maps of polarized synchrotron 
emission of the SNR. 

2.2 Geometry 

Figure [T] shows a 3-D sketch of the simulated SNR and 
its projection along the line of sight. To simulate the im- 
ages of the SNR shell we assumed that a quarter part of a 
spherical forward shock has a relativistic electron distribu- 
tion N(z,j,t) that does not depend on the azimuthal and 
polar angles, but is inhomogeneous in the radial direction 
with a strong peak (of width A s ) at the shock position at 
r — -Rsnr- The line of sight is along the z axis. The Stokes 
parameters I, Q, U, V for incoherent photons are additive, so 
we can integrate these over the line of sight weighted with 
the distribution function of emitting particles N(z, 7, t) to 
get the intensity 

S(R_L,t,v) = J dzd<yN(z,~f,t') J(r,7,t'). (1) 

To collect the photons reaching the observer at the same mo- 
ment t, we performed an integration over the source depth 
using the retarded time t' = t — |r — Rj_|/c as argument in 
B(r, t') and N(r, E,t'). The integration grid has a cell size 
smaller than L m i n . The result is a surface density of Stokes 
parameters of radiation from the volume along the line of 
site. The fourth Stokes parameter V is zero in the case of an 
isotropic electron velocity distribution. In order to achieve 
a few percent accuracy we integrated over 8000 grid points 
along the line of sight. The number of pixels in the sky pro- 
jection is 100 x 200. The degree o f polar ization was derived 
following iGinzburg fc Svrovatskiil (| 19651 ). 

Below we present synchrotron images simulated with a 
steady model distribution of electrons accelerated by a plane 
shock of velocity 2,000 km s _1 propagating in a fully ion- 
ized plasma of number density 0.03 cm -3 . The kinetic model 



used to s imulate the electro n distribution was described in 
detail bv lBvkov et~aH (|2000l ). The magnetic field in the far- 
upstream region was fixed at 3 fiG and it was assumed that 
the magnetic-field amplification produces a random field of 
sj (B 2 ) — 3 x 10~ 5 G in the shock vicinity. In the ran- 
dom magnetic field simulations we used a wavenumber range 
fcmin < k < fc max , where L min = 27r/fc max = 2 x 10~ 4 7r • D 
for 8 = 1.0 and L min = 27r/fc max = 2 x 10~ 3 tt • D for 8 = 2.0, 
with L max = 27r/fc m i n = 0.2-7T • D for both 8 values. Here D 
is the size of a unit cubic box, as shown in Figure [1] The 
random magnetic field in the simulated SNR shell quarter 
was divided into 8 such boxes. This number of boxes was 
chosen to achieve the required accuracy of the integration 
of the random field along the line of sight. The field in the 
boxes was simulated as a function of global SNR coordinates 
as it is shown in Figure [1] (i.e., not just locally in each box). 
Note that the field was actually simulated in a region larger 
than the SNR shell and that the box sizes are larger than 
the sizes of the random filamentary structures that appear. 



3 SIMULATED POLARIZATION MAPS OF 
THE X-RAY SYNCHROTRON EMISSION 

Figures [2]|4] show examples of the resulting maps at differ- 
ent X-ray energies. The left panels show the synchrotron 
intensity, the right panels the polarization degree, and the 
central panels the product of the two. The latter is a mea- 
sure of the polarized flux and is meant to illustrate that 
peaks in the polarization-degree map do not necessarily cor- 
respond to peak intensities. The images clearly demonstrate 
1) the presence of detailed structures - clumps and filaments 
- produced by the stochastic field topology (for details see 
iBvkov et ai]|2008ft and 2) that some of these structures emit 
highly polarized emission (> 30%) at energies of 5 keV and 
above. 

Figure [5] shows the 5 keV map (as in Figure |3J, but for 
a steeper spectrum of magnetic fluctuations (6 = 2.0 rather 
than 1.0). There is a distinct difference, indicating that for 
steeper spectra the size of the polarized structures is larger 
and the degree of polarization of these structures is higher 
(about 50% for 8 = 2.0 and *JWj = 3 x 10" 5 G ). 

Irrespective of the precise value of 8, it is clear from Fig- 
ures [2][5] that the degree of polarization is higher at higher 
X-ray energies. The physical reason is best illustrated in 
case of a power-law electron distribution (with spectral index 
r). Namely, the degree of polarization II and the local syn- 
chrotron cmissivity I(r, t, v) have the following dependencies 
on V: II w (r + l)/(r + 7/3) (i.e. the degree of of polarization 
IT is increasing with Y) and /(r, t, v) tx 5 1 / 2 ( r+1 ) (i.e. the lo- 
cal emis sivity is relatively very high for large B and large F) 
(see e .g., IGinzburg fc Svrovatskiil 19651 ; iRvbicki fc Lightmanl 
Il979l ). In the high-energy cut-off regime the electron spec- 
trum is typically exponential, but the effective index V is 
large indeed and the value increases for electrons emitting 
at higher frequency v. This explains the increase of the po- 
larization degree with v. In addition, the dependency of the 
emissivity I on B and V can lead to a highly polarized bright 
feature that stands out in the map for even a single strong lo- 
cal field maximum. In lower-energy maps, for which V on av- 
erage is smaller (i.e. well below the cut-off regime), high po- 
larization of a single maximum can be smoothed or washed 
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Figure 2. Simulated maps of polarized synchrotron emission in 
a random magnetic field at 0.5 keV. Intensity, v 2 ■ I(R±,t, u), is 
shown with a linear color scale in the left panel. The central panel 
shows the product of intensity and polarization degree. The right 
panel shows the degree of polarization indicated by the colorbar. 
The stochastic magnetic field sample has »/ (B 2 ) = 3 X 10 -5 G 
and spectral index 6 = 1.0. 




0.1 0.15 



0.25 0.3 



Figure 3. The same maps as in Figure[2] but at 5.0 keV. 

out by contributions from a number of weaker field maxima 
integrated over the line of sight. This effect can be seen in 
Figures r2l4l The high-energy maps are 'twinkling' because of 
the finite life-time of the magnetic-field amplifications. The 
timescale for variations in the polarization (and the energy 
dependence) is similar to that of the time variability of the 
intensity maps (studied in §4.1 of lBvkov et al.ll200&f ). 

Figure |S] illustrates the dependency of the polarization 
degree on the resolution of the simulated maps (9", 18", 
and 36") at 5 keV. In Figures [7] and [5] the polarization maps 
are presented at 20 keV for larger pixel sizes of 3' and 7.5' 
(close to the INTEGRAL ISGRI pixel size) and 8 = 1 and 
2. Comparison of Figures[7]and[5]shows again the strong de- 
pendence on the spectral index S of the stochastic magnetic 
field. 
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Figure 4. The same maps as in Figure[2] but at 50.0 keV. 




0.05 0.1 0.1 5 0.2 0.25 0.3 0.35 



.45 0.5 0.55 



Figure 5. The 5.0 keV synchrotron map for a different magnetic 
field spectrum in the shell. The stochastic magnetic field sample 
has \J(B 2 } = 3 X 10~ 5 G and 5 = 2.0. 




Figure 6. The simulated 5.0 keV synchrotron polarization maps 
with different pixel sizes. The left one has a pixel size of 9", the 
central of 18", and the right of 36". The yellow frame of 2.6' X 2. 6' 
indicates the field of view of the XPOL polarimeter (see text). The 
stochastic magnetic field sample has y/ (B 2 ) = 3 X 10 — 5 G and 
5 = 1.0. The simulated SNR shell has the radius of about 0.4°. 



4 OBSERVATIONAL PERSPECTIVE 



This work was stimulated by the fact that a number of X-ray 
polarimeter instruments is being considered currently. The 
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Figure 7. The simulated 20.0 keV synchrotron polarization 
maps with different pixel sizes. The left one has a pixel size 
of 3' and the right of 7.5' . The yellow line indicates the for- 
ward shock position. The stochastic magnetic field sample has 
^/{B 2 ) = 3 X 10~ 5 G and S = 1.0. The field was simulated in a 
box larger than the SNR shell, but the regions well outside the 
forward shock arc dim as it is clearly seen in the left panels in 
Figures [2]- [5] 
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Figure 8. The same maps as in Figure [7] but for magnetic 
turbulence with S = 2.0. 



Im aging X-ray Polarimet ry Explorer IXPE was proposed 
bv lWeisskopf et al.l (|2008h as a dedicated X-ray-polarimetry 
observatory to measure the X-ray line ar polarization as a 
function of energy, time, and position. iLegere et al.l (|2005l ) 
is developing a Compton polarimeter to measure polariza- 
tion of hard X-rays in the 50-300 keV energy range. A 
ba lloon-borne hard X-Ray p olarimeter HX-POL is proposed 
bv lKrawczvnski et all (|2008ft . A Hard X-ray Telescope HET 
aboa rd the Energet ic X-ray Imaging Survey Telescope EX- 
IST |Grindlavll200l ), with a wide field of view for the coded 
aperture imaging, is being designed to study the polariza- 
tion at high energies and its temporal evolution. Polariza- 
tion detection of X-ray sources as faint as 1 milliCrab is an 
aim of the Gravity and Extreme Magnetism SMEX (GEMS) 
mission that u ses foil mirrors and Time Projection C ham- 
ber detectors (|Swank et al.l 120081 : [jahoda et al.1 120081 ). The 
missions listed above have good perspectives in this field 
of research. We address here the p otentials of the X POL 
polarimeter as proposed for XEUS (|Costa et alj|2008h - al- 



though evolved into International X-ray Observatory (IXO) 
in the mean time - to illustrate the observational possibili- 
ties for synchrotron X-ray studies of SNR RXJ1713.72-3946 
as a generic example. 

RXJ1713. 72-3946 has an extended shell of about one 
degree angular diameter. A field of view of 1. 5' x 1.5' and 
an an gular resolution of 5" was proposed by ICosta et al.l 
{2(xjl), if XPOL was part of the XEUS mission. A polarime- 
ter like XPOL aboard the IXO mission Q will have a some- 
what larger field of view. We illustrate that case in Figure [S] 
where the field of view of 2.6' x 2.6' is shown as a yellow box 
in the polarization maps simulated with pixel sizes of 9", of 
18", and of 36" to illustrate the angular resolution effect. 

In the case of an extended source like RXJ1713.72- 
3946 first of all a wide field map of the source region 
is needed to specify the XPOL pointing. Using Chandra 
archive data, we estimate the 2-10 keV flux from a 2.6' x 
2.6' region in the shell of RXJ1713. 72-3946 to be about 
4.5x 10 -12 erg cm~ 2 s~ . From the minimum detectable 
polarization as a function of o bserving time as presented 
in Fig. 6 of ICosta et al.l (|2008h . corrected for the reduced 
effective area of IXO 1 , we estimated that a meaningful po- 
larization map can be constructed with XPOL within an 
exposure time of 100 ks. The polarization map of 2.6' x 2.6' 
should likely reveal detailed highly-polarized structures with 
typical sca les of about 10" as it w as discovered in Chandra 
images by lUchivama et al.l (|2007l ). The degree of polariza- 
tion will increase with increasing X-ray energy as predicted 
from the modeling in this paper. The polarization will be 
time variable on a few year time s cale (dependi n g on the 
photon energy) as it was found by iBvkov et al.l ( 20081 ) in 
simulated intensity maps. Mapping of the whole extended 
shell of RXJ1713. 72-3946 in polarized X-rays would be un- 
realistic with that set up, therefore the target must be first 
identified with a wide field X-ray imager. 

Another object of great interest is Cas A, that is of an 
angular size comparable with the field of view of XPOL. 
We estimate that some thin peripheral polarized X-ray fil- 
aments of some ten arcsecond scale can be studied with 
XPOL, also with an exposure of about 100 ks. In the DSA 
model the scale size L max of the magnetic fluctuations re- 
sponsible for the twinkling polarized structures is expected 
to be connected to gyroradii of accelerated protons at max- 
imal energies. These can be roughly estimated to be about 
3 x 10 17 ■ B~ G ■ EiooTeV cm. Therefore, their angular sizes 
are expected to be above a few arcseconds for SNRs within 
a few kiloparsec distance. 



5 CONCLUSIONS 

We have studied the polarization of X-ray synchrotron emis- 
sion from SNRs addressing the significant effect of magnetic- 
field fluctuations on synchrotron emission in X-rays. Such 
magnetic fluctuations form a natural starting point because 
they must be present if diffusive shock acceleration is indeed 
the basic me c hanis m for accelerating particles in SNRs. Like 
IBvkov et al.l (|2008h we simulated random magnetics field to 
construct synchrotron emission maps, given a smooth and 



http: / / ixo.gsfc. nasa. gov / science / performanceRcquircments .html 
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steady distribution of electrons, but now with special at- 
tention to the polarization of the resulting emission. The 
simulated random magnetic fields show non-steady local- 
ized structures with exceptionally high magnetic-held am- 
plitudes. These magnetic-field concentrations dominate the 
synchrotron emission - integrated along the line of sight - 
from energetic >TeV electrons, i.e. in the cut-off regime. In 
terms of a power-law electron spectrum with spectral index 
r, this can be understood since the synchrotron emissivity 
7(r, t, v) is proportional to 5 1 / 2 ( r + 1 ) (i.e. the local emis- 
sivity is relatively very high for large B and large F). The 
power-law approximation is only useful over a narrow elec- 
tron energy range in the cut-off regime, where the effective 
spectral index F is increasing with the electron energy. 

Starting from the simulated magnetic-field simulations, 
we have constructed maps of polarized X-ray emission of 
SNR shells. These are highly clumpy with high polarizations 
up to 50%. This characteristic of high polarization again ap- 
plies to energetic >TeV electrons in the cut-off regime. In 
terms again of a power-law electron spectrum with spectral 
index F, this can be understood since the degree of polariza- 
tion is given by II w (r + l)/(r + 7/3) (i.e. II is increasing 
with r). 

The distinct characteristic of the modeled synchrotron 
emission is its strong intermittency, directly resulting from 
the exceptionally high magnetic-field amplifications ran- 
domly occurring as shown in the simulations. Also charac- 
teristic is the increase of the polarization degree with X-ray 
energy addressed in |3] Since this "twinkling" polarized X- 
ray emission of SNRs depends strongly on the magnetic-field 
fluctuation spectra, it provides a potentially sensitive diag- 
nostic tool. 

The intermittent appearance of the polarized X-ray 
emission maps of young SNR shells can be studied in detail 
observationally with imagers of a few arcsecond resolution, 
though even arcmin resolution images can provide impor- 
tant information as it is illustrated in Figures 161 7181 The 
polarized emission clumps of arcsecond scales are time vari- 
able on a year or longer (depending on the observed photon 
energy, magnetic field amplification factor and the plasma 
density in the shell) allowing for rather long exposures even 
in the hard X-ray energy band. Hard X-ray observations 
in the spectral cut-off regime are the most informative to 
study the magnetic fluctuation spectra and the acceleration 
mechanisms of ultra-relativistic particles. 

Altogether, the modeled appearance and its time vari- 
ability - on a timescale of typically a year - resembles closely 
what is observed already in X-ray images of some young 
supernova remnants. Observing the predicted high polar- 
ization in clumps and filaments, however, should proba- 
bly await future instruments that are currently being con- 
sidered. Such observations will provide unique information 
on magnetic fields and high-energy particle acceleration in 
SNRs. 
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